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Abstract 

In this paper we study a class of stationary states for reaction-diffusion systems 
of > 3 densities having disjoint supports. For a class of segregation states governed 
by a variational principle we prove existence and provide conditions for uniqueness. 
Some qualitative properties and the local regularity both of the densities and of their 
free boundaries are established in the more general context of a functional class char- 
acterized by differential inequalities. 

1 Introduction 

The occurrence of nontrivial steady states (pattern formation) for reaction-diffusion sys- 
tems has been widely studied in the literature. Of particular interest is the existence of 
spatially inhomogeneous solutions for competition models of Lotka-Volterra type. This 
study has been carried out mainly in the case of two competing species, see e.g. [^ [inilT^ 
Em E21 ESI EHl 123 EH] ; in recent years also the case of three competing densities, which is 
far more complex, has become object of an extensive investigation jl3[ll4lOHl25j . In most 
cases, the pattern formation is driven by the presence of different diffusion rates when the 
coefficients of intra-specific and inter-specific competitions are suitable related. A remark- 
able limit case of pattern formation yields to the segregation of competing species, that is, 
configurations where different densities have disjoint habitats; see [TTl [T^ [TBI [T71 [T5ll^ . 
Object of the present paper is to study a class of possible segregation states, involving an 
arbitrary number of competing densities, which are governed by a minimization principle 
rather than competition-diffusion. Roughly speaking, we are going to deal with stationary 
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configurations of A; > 2 densities tliat interact only through the boundaries of their nodal 
sets; the minimization involves the sum of the internal energies, with the constraint of 
being segregated states. In other words, the supports of the densities have to satisfy a 
suitable optimal partition problem in M.^ . Precisely, let Q be a bounded open subset of 
{N > 2) and let us call segregated state a /c-uple U = {ui, . . . ,Uk) S {H^{0,))^ where 

Ui{x) ■ Uj{x) =0 i j, a.e. x E il. 

We define the internal energy of U as 

J{U)= J2 y^(^ld^{x)\Vui{x)\^ -Fi{x,Ui{x))^dx 

where the Fj's satisfy suitable assumptions (see (Al), (A2) below). Our first goal is to 
minimize J among a class of segregated states subject to some boundary and positivity 
conditions; next we shall develop a regularity and a free boundary theory for minimizers. 
In performing the second goal the main tools will come from recent results that Caffarelli, 
Jerison and Koenig (see also [HISl El El and references therein) have obtained in the 
study of free boundaries in other contexts. Recently, free boundary problems have been 
studied in connection with the asymptotic behaviour of some models of population dy- 
namics with diffusion as in Dancer and Du (T^) and Dancer, Hilhorst et al in ( JZ|)) in the 
case of two competing species. In a forthcoming paper we shall show how our variational 
problem appears as a limiting problem for some classes of competition-diffusion problems. 

Our first result establishes existence for this problem; then we discuss the uniqueness of 
the solution. Surprisingly enough, the minimizer can be proven to be unique for a large 
class of Lagrangians. A remarkable fact is that solutions to this variational problem satisfy 
extremality conditions in the form of differential inequalities of special type. Precisely they 
belong to the functional class 

> 

Ui > 0, Ui ■ Uj = i( i ^ j 
-Aui < fi{x,Ui), -Aui > f{x,Ui) 

where Ui = Ui - Y^h^i Uh and f{x,Ui) = J2j fji^, Ui)x{u,>o}- 

A further reason of interest in the class S is that it contains also the asymptotic limits 
of the solutions of a large class of competition-diffusion systems when the inter-specific 
competition terms tend to infinity. This will be the object of a forthcoming paper; a link 
between some variational problems and competing species systems has been traced by the 
authors in [HI EI • 

An important part of the paper is devoted to study the qualitative properties exhibited 
by the segregated states belonging to the class S. In particular we shall establish the 
local lipschitz continuity both of U and its nodal set; to this aim we will take advantage 
of some monotonicity formulae as in [21 0- Then, for the dimension = 2, we develop 
further our investigation. Our main result is that, near a zero point, U and its null set 
exhibit the same qualitative behavior of harmonic functions and their nodal sets ([ni21j). 
In particular we prove that the set of double points (i.e. points where two densities meet) 



{ui, 



,Uh, 
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is the union of a finite number of regular arcs meeting at a finite number of multiple points 
(i.e. points where more than two densities meet). We emphasize that, at a multiple point 
the densities share the angle in equal parts and moreover an asymptotic expansion for U 
is available. 

The plan of the paper is the following: in Section ^ we introduce the basic assumptions 
and formulate the variational problem; the existence of a minimizer is proven in Section 
13 Section |1] deals with the uniqueness of the solution; finally in Section [21 the extremality 
conditions are established. In Sections El and [3 we introduce S and a wider functional 
class S* ; next in Section |H1 we prove the local lipschitz continuity in S* and the global 
regularity in S. In Section (HI we establish some qualitative properties of the elements of S 
in dimension N = 2. 

2 Assumptions and notation 

Let TV > 2; let 17 C be a connected, open bounded domain with regular boundary dQ. 
Let A: > 2 be a fixed integer. Throughout all the paper we will make the following set of 
assumptions (for every i = 1, . . . , k): 

• 4>i €z H^^^{d^l), (f)i > 0, (pi ■ (f>j = for i ^ j, almost everywhere on d^l; sometimes 
such a boundary datum will by called admissible 

• die W^'°°{Q), di > onH 

• fi{x, s) ■MxR+ such that: 

(Al) fi(x,s) is Lipschitz continuous in s, uniformly in x and fi{x,0) = 
(A2) there exists bi e L°°{U) such that both 

\fi{x, s)\ < bi{x)s Vx € ri, s > s » 1 

and 

J (^df{x)\Vw{x)\^ -bi{x)w'^{x)^dx>0 ^weH^in). 

Every (j)i will be the boundary trace of a non negative density Ui G H^{Q). Moreover, 
associated to each density, we consider its diffusion coefficient di and its internal potential 
Fi{x,s) := fi{x,u)du. 

We are concerned with the following variational problem. 
Problem 2.1 Let 

U = |(iii, ■■■ ,Uk) G {H^{^))'' : UilsQ = c/)i,Ui > Oy i = I, - ■ ■ ,k; uj ■ Ui = 0,i ^ j a.e. on Q 
Find the minimum of the functional 

'^(f^)= E [J^(^ldi{x)\Vu,{x)\'-F,ix,u,ix))yx'j (1) 

where U eU. 
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Remark 2.1 We notice that the case with variable diffusions can be reduced, by a suitable 
change of the variables Uj's, to the case when di = 1 for every i. Indeed, let Ui = Vi/di 
(this can be done because the dj's are strictly positive). After integrating by parts we 
obtain 

= E -Fi{x,Vi)^dx-C^ = J{V)-C. 

Exploiting this identity, the reader can easily check that the validity of the assumptions 
for the /j's and for the (/>'s implies the same for the new data /j's and (/)i's. Hence, in what 
follows, we will choose di = l for every i, unless otherwise specified (namely in the results 
of Section EJ. 

Remark 2.2 By our definition, the functions /j's are defined only for non negative values 
of s (recall that our densities ttj's are assumed non negative); thus we can arbitrarily define 
such functions on the negative semiaxis. For the sake of convenience, when s < we will 
let fi{x,s) := —fi{x,—s). This extension preserves the continuity, thanks to assumption 
(Al). In the same way, each Fi is extended as an even function. 

Notation In the following, when not needed, we shall omit the dependence on the variable x. 
We use the standard notation g^{x) = maxj;^Q{g{x),0) and g^{x) = maxx^^i{—g{x),0). Given a 
fc-uple (ui, . . . , Ufc) we introduce the "hat" operation as 

:= Ml - Uj. 

Furthermore, with some abuse of notation, we shall use a capital letter to identify both a /c-uple 
and the sum of its k components (e.g. U = (mi, . . . ,Uk) and U = X^iLi )■ With the notation 
iui,n) we shall denote the i-th component of a sequence of fc-uples (C/„). The symbol xa will 
denote the characteristic function of the set A. 



3 Existence of the minimum and continuous dependence 

Our first goal is to prove the existence of at least one minimizer of Problem 12. II Next we 
shall prove continuous dependence of the minimizers (that do not need to be unique) with 
respect to the data. This continuity property will be exploited in the analysis of the local 
properties of the solutions, and prescisey when performing the blow-up argument. 

To this aim, we start observing that our assumptions on fi imply 

\F,{x,s)\<^^s^ + C\s\ Vxer?, VsgM (2) 

for every i. On the other hand, by standard eigenvalues theory, assumptions (A2) implies 
that the quadratic form there is an equivalent norm on //^(il), that is, there exists e > 
such that 

/ {\Vw\'^ - bi{x)w^)dx > e [ \Vw\^dx, 
Jn Jn 
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for every w € Hq{^) and for every i. As a consequence we have the following result: 

Theorem 3.1 Under the assumptions o f Section Wi Problem 12 . l1 has at least one solution. 

Proof: applying Q we easily obtain that the if^-continuous functional is coercive; 
indeed, for every Ui G H^{^), we have 

(^|Vni|2 - Fi{u,)^ dx> jj^ (iVuip - biufj dx- Cmdx 

> J / |Vn,j2-c(l^) 
2. Jn 

for some constant c{Q). Let us take a minimizing sequence (ui,n) in U of disjoint sup- 
port functions. This sequence being if^-bounded by the above inequality, there exists a 
subsequence weakly convergent to Uj in and, by compact injection, in the L^-strong 
topology; taking possibly a new subsequence, we infer almost everywhere convergence in 
n of every Ui^n to Ui and the limit functions have obviously disjoint supports. The weak 
lower semicontinuity ensures that the weak limit is in fact a minimizer. ■ 



Theorem 3.2 Let Un = (wi,n, ■ • ■ ,Uk,n) be a sequence of solutions to Problem l2.1L respec- 
tively with admissible data {(pi^n, • • • , 4'k,n), such that 

(t>i,n^(t>i inH'/\dn) 

and potentials 

Fi,n{x,s) — >Fi{x,s) inC^inxR) 
for every i = 1, ■ ■ ■ ,k. Then, 

Ui,n — ^ Ui inH^{Q) 

for every i = 1, ■ ■ ■ , k and U = {ui, • • • , Uk) solves Problem 12.11 with data {(pi, ■ ■ ■ , cpk) cmd 
potentials {Fi, ■ ■ ■ , Fk). 



Proof: first, we observe that {(pi, ■ ■ ■ , (j)k) is an admissible datum, i.e. the <j)j^s are 
nonnegative and have disjoint supports by the strong convergence of {(pi,n-, - ■ ■ ,(pk,n) in 
H^l^{d^). 

We denote by U* = {u\, ■ ■ ■ ,u*^) a solution to Problem 12.11 with data {(pi, ■ • • ,(pk), and 
{Fi, ■ ■ ■ , Fk). Consider the minimum levels 

c„ = J{Un) and c* = J{U*). 

Observe that the convergence of the boundary traces (pi^nS and of the Fi^^s, ensures a 
bound on the sequence c^. The coercivity of J then yields to a bound on the sequence 
ll^j,n||_ffi(n); therefore we can assume, up to a subsequence, that 

Cji > c 

Ui^n — ^ u-i weakly in H^{n). (3) 



5 



Furthermore, as a consequence of the compact injection ^ L^, we have Uj ■ uj = 0, 
whenever i ^ j. Moreover, by the weak continuity of the trace operator, we obtain 

Ui\dn = 

The lower weak semicontinuity of the norm imphes 

- V / |VMip<hminf- V / |Vni„P 

and we also have 

/ Fi{x, Ui{x)) = Ihn^ / Fi^ri{x,Ui^n{x)). 



We observe that the level 

c 



j=l---fc 

is not necessary a minimum level but satisfies, by the discussion above, the inequalities 

Co > c> c*. 

We wish to prove that cq = c*. Suppose, by contradiction, that c* < cq. Consider the 
harmonic extensions (still denoted with the same symbols) on 0, of the (pi^nS and of the 
cpi^s and introduce ipi^n = 4'i,n — <t>i- Then, by construction 

^,,„ inH\^) (4) 



We define 

Wi,n 
Vi n 



We observe that Wi^nldn = moreover, since u* > and by @, 

^ u*mH\n). (5) 

Moreover, since tfj^n > 0, it is immediate to see that Vi^n ■ Vj^n = Hi ^ j- Hence it follows 
from the definition of Cn that 



on the other hand © implies 

_E ji/jV.,,„l^-/^F.,„K„)}^c- 

that implies that c* > cq. 

Finally, from the equality of the minima cq = c = c* , we also deduce the strong convergence 
of the Ui^nS in H^{i}) and the thesis follows. ■ 
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4 Uniqueness 

In general we can not expect the minimizer of Problem 12.11 to be unique. A simple 
counterexample can be constructed in the following way: 

Example 4.1 Let k = 2, (pi = (p2 = 0, fi{x,s) = f2{x,s) = min{\s,s^/^), where 
A > Ai(r2), the first eigenvalue of the Laplace operator with Dirichlet boundary condi- 
tion. The reader can easily check that the assumptions of Theorem \3.1\ are fulfilled and 
that the infimum is smaller than zero. Moreover the infimum of the associated functional 
is achieved by a function of the form {u, 0), or (0, u), where u ^ 0; indeed, any minimizer 
(ni,M2) can be replaced by, say, {ui + U2,0) keeping the same level of the functional. Hence 
the associated variational problem does not have uniqueness of solutions. 

A major obstruction to uniqueness is the lack of convexity that may occur both in the 
Lagrangian (as the above example illustrates) and in the constraint; nevertheless, the 
following result shows that the full convexity of the Lagrangian is sufficient to prove 
uniqueness of the minimizer, provided the diffusions do not depend on i. 

Theorem 4.1 Under the assumptions of Section\^ assume moreover that 

di = dj ,yi,j (6) 

Fi{x,s) is concave in the variable s, for all x ^Q. (7) 
Then Problem 12.11 has an unique minimizer. 



Proof: let c = inf{ J(C/) : U G lA}. Arguing by contradiction, we consider two minimizers 
U = {ui, • • • , Uk) and V = (vi, • • • , Vk) achieving c, with Uj ^ Vi for some i. For every i 
and A S [0, 1] , we define 

Ui = Ui-^Uh (8) 
Vi = Vi-^Vh (9) 

hj^i 

wl^^ = [Xui + {l-X)v,]+. 

Our goal is to show that J(w|'*'^) < XJ{ui) + (1 — X)J{vi), for every A € (0,1). It is 
worthwhile noticing that this property can be seen as a convexity type property combined 
with a special type of projection on the constraint U. 

To begin with, we have to show that the w^^^^s satisfy the constraint. We first notice that 
u!i^^ > and WiIqq = (pi. Furthermore we have 

w^^^vUj^'^ = a.e. on Q , for j / i. 

Indeed, assume wl'^\x) > 0; this means that 

Xui{x) + (1 — X)vi{x) > ^ Xuh{x) + (1 — X)vh{x) > Xuj{x) + (1 — X)vj{x) Vj ^ i. 
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Therefore we have, when j ^ i, 

Xuj{x) + (1 — X)vj{x) < \ui{x) + (1 — X)vi{x) < ^ Xuh{x) + (1 — X)vh{x) , 

and hence Wj^\x) < 0. 
Let us denote the supports 

rf ) = {-f ^ > 0} ; 
recalhng that Fi{x,0) = we note that 

i=l i=l •'^t 

In view of ©, using the convexity of the quadratic part of the functional, and the defini- 
tions ((HI),© and keeping in mind that both the {u-i)' s and {v-i)' s have disjoint supports, 
we obtain that, for every A G (0, 1), 

k 

^ /^(.) (id2|Vu;W|2)dx < Eii /rw A(id2|v2,|2) + (i _ X){\<p\Vv,\'')dx 

1 = 1 

< Eii In A(id'|Vn.|2) + (1 - X)iy'\Vv.\')dx 

Now we turn to the potential integral. By assumption ((T)) and the evenness of the potentials 
Fj's, the inequality 

-Fi{x,Xui{x) + (1 - X)d,{x)) < -XFi{x,Ui{x)) - (1 - X)Fi{x,Vi{x)) 

holds whenever x G {n,, > 0} fl {vi > 0}. Hence, assume that x € {ui > 0} fl {vj > 0}, for 
some j 7^ i; let us fix A and let x such that Xui{x) — (1 — X)vj{x) > (the symmetric case 
is obtained by parity), and we study the term —Fi{x, Xui{x) — (1 — X)vj{x)). Introducing 
the auxiliary function ^'(A) = —Fi{x,Xui — (1 — X)vj) + XFi{x,Ui) + (1 — X)Fj{vj), let 
A G (0, 1) such that Xm - (1 - X)vj = 0. It is easy to see that ^'(1) = 0, and ^'"(A) > 0. 
Moreover, ^'(A) < 0, since Fj(x,0) = and,by Q, Fi < 0. Hence, by convexity, we infer 
that ^'(A) < 0, for every A S (A, 1) and therefore 

-Fi{x, {Xui + (1 - A)%)+) < -XFi{x,Ui) - (1 - X)Fj{x,Vj) 

holds when x € {ui > 0} fl {vj > 0}. 

Finally, gathering together all these inequalities, for every fixed A we obtain 
J( Y: wi'^) < A J( Y: u.) + (1 - A) J( ^ t;.) < c , 

i=l i=l i=l 

a contradiction. ■ 
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The requirement of Theorem l4.1l that the Lagrangians are convex in ah variables may seem 
very restrictive. The following result makes a different assumption, still sufficient for the 
uniqueness, that may be more useful in the applications, for it is always satisfied (given 
the dj's and -Fj's) provided the domain Q is small enough. 



Corollary 4.1 Let di e W^'°° be given. Assume that Fi are of class C in the variable s 
and let us denote 

bi{x) = sup-^{x, s) (10) 
Assume that there exists a positive function d such that, for every i = 1, . . . ,k, 

Ad+(^-^]d>0 ini}. (11) 



Then, Problem 12.11 has a unique minimizer. 



Proof: the following identities hold: 



dl(x)\Vui\^ 



d^(x) 



d^ix) 



V 



V 



diUi 



d 



+ div ( f^'V ( ^ 



d J 

d 



an 



ddiufV ( ^ I • 



Therefore, up to a constant (recall that our boundary data are prescribed), by the change 
of variables Vi = diUi/d we can transform the initial Lagrangians into ones of the following 
form: 

d\x)\Vv,\^ + ^div (d^V (%\] vf - Fi{x,dvi/di). 



di \ \ d , , 

Recalling the definition of the 5j's of (|1()|) . one easily checks that these Lagrangians are 
convex in the variable Vi provided the following inequality holds for every index i 

d^ 



di V V a 
that is easily seen to be equivalent to (fTT|) . 



,b, > 



5 Extremality conditions 

The goal of this section is to prove that the minimizers of the variational problem Problem 
12.11 satisfv a suitable set of differential inequalities. To start with, let {ui, . . . ,Uk) be a 
minimizer of problem 12. H we define f{x,Ui) 



f{x, Ui) = ^ fj{x, Ui)X{u^ 



>o} 



fi{Ui 



if x G {ui > 0} 



-fj{x,Uj) X e {uj > 0}, j ^ i. 



(12) 



9 



Note that this definition is consistent with that of (jH)), for the functions fi are extended 
by oddness. Our main result is the fohowing: 

Theorem 5.1 Let U be a solution to Problem I2.1L Then, for every i, we have, in distri- 
butional sense, 

(i) -Aui < fi{x,Ui) 
(a) -Aui > f{x,Ui) . 

Proof: 

(i) We argue by contradiction. Then, there exists at least one index j such that the 
claim does not hold; that is, there is < € C'^{Q.) such that 

j VujVcj) - fj{x, Uj)(j) > 0. 

For < t < 1 we define a new test function V = (fi, . . . , f^) as follows: 

{Ui if i / j 

{ui - t(t))^ if i = j. 

We claim that V lowers the value of indeed we have 

J{V) - J{U) = IJ (|V(n,- - - |Vn,f ) - ^F,(x, {uj - t4>)+) - F,{x,u,) 

- Jul {I'^Uj - t(t>\^ - \VUj\^) +t Jjj{x,Uj)(f> + 0{t) 

< t [ {-VujVcj) + fAx.uM) + o(t). 

Choosing t sufficiently small, we obtain 

J{V) - J{U) < 0, 

a contradiction. 

(ii) Let j and < (/) G C^{^) such that 

/ SJujSJcj) — f{x, Uj)(f> dx < 0. 

Again, we show that the value of the functional can be lessen by replacing U with an 
appropriate new test function V. To this aim we consider the positive and negative 
parts of Uj + tcj) and we notice that, obviously, 

{{uj + 14>)- > 0} C {{uj)- > 0} = Ui^j{ui > 0} . 
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Let us define V = {vi, . . . ,Vk) m the following way: 

{uj + tc/))^ , if i = j 

We compute, using the definition ((T^ . 



J{V)-J{U) = Yl /^^(IV^.P-IV^X, 



- {Fi{x,Vi) - Fi{x,Ui)) 

\ (|V% + t^f - |V%|2) - (Fj(x, (uj + t^)+) - Fj{x,uj)) 
J2 {Pi {x, {uj + tcpy Xk>o}) - Fi{x,Ui)^ 

t VUjV<j)- fj{x,Uj)x{u,>o}4' + y2fi(^^Ui)x{u,>o}4' + o(t) 



= t / VujVcf)- / fj{x,Uj)(l) + o{t) . 
Jn Jn 

For t small enough we find J{V) < J(U), a contradiction. 



6 The class S and its basic properties 

Let U be the set of admissible A:-uples as defined in Problem 12.11 and let fi be given 
satisfying (^1), (^2); we introduce the following functional class: 

Definition 6.1 

5 = <^ (m, ■■■ ,Uk) eU : -Aui < fi{x,Ui) , -Aui > J{x,Ui),\/i = 1,... ,k 



By virtue of Theorem l5. li the class S is the natural framework where to develop our theory 
of regularity and free boundary of minimizers to Problem 12.11 As already noticed in the 
introduction, the class S is of independent interest, for it contains the asymptotic limits 
of highly competing diffusion systems. 

Let us start with the following definitions: 

Definition 6.2 The multiplicity of a point x £ Q is 

m{x) = jj {i : meas {{ui > 0} fl i?(x, r)) > Vr > 0} . 
We shall denote by 

Zh{U) = {x e O : m{x) > h} 
the set of points of multiplicity greater or equal than /i G N. 
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The following properties are straightforward consequences of the definition of 5 joint to the 
locally Lipschitz continuity of the /j's that implies the validity of the Maximum Principle 
for elliptic equations. 

Proposition 6.1 Let x € Jl.' 

(a) If m(x) = 0, then there is r > such that Ui = on B{x,r), for every i. 

(h) If m{x) = 1, then there are i and r > such that Uj > and —Aui = fi{x,Ui) on 
B{x, r). 

(c) If m{x) = 2, then are i,j and r > such that = for k ^ i,j and —A{ui — uj) = 
fij{x,Ui -Uj) on B{x,r), where fi,j{x,s) = - fj{x,s~). 

Proof: part (a) follows directly by the definition of multiplicity; assume m(x) = 1 and let 
r > be such that meas{{uj > 0} H B{x, r)) = for every j except i. Then, by definition 
of 5, Ui satisfies the equation — A-Uj = fi{x,Ui) on B{x,r). We can write 

-Aui = ai[x)Ui = r~^'^i 

Ui{x) 

where Oj E L°° by the Lipschitz continuity of the /j. Then, since Uj > 0, we infer from 
the strong maximum principle that > on B{x,r). The second statement follows 
immediately from the definition. ■ 



Remark 6.1 

- We can not exclude, at this stage, the occurrence of points of multiplicity zero, although 
this possibility will be ruled out at the end of Section 19.21 at least in two dimensions, 
under a weak non degeneracy assumption. Note that d{x £ O : m(x) = 0} C ^3 U di}. 

- The second point of Proposition I6.1l savs in particular that, if m{x) = 2 

lim \/ui{y) = - lim Vuj{y) . 

t/e{uj>o} !/e{iij>o} 

If, by the way, the above limit is not zero, it follows that the set {x : m(x) = 2} is locally 
a manifold of dimension — 1. Of course the above equality has to be changed, in 
presence of variable diffusions, into 

1™ di{x)Vui{y) = - lim dj{x)Vuj{y) . 

A major goal in the subsequent analysis will concern the geometrical properties of the 
supports of the densities Ui and their common boundaries. As a first consequence of the 
Maximum Principle we can give a criterium for the connectedness of the supports: 

Proposition 6.2 Assume, for some i, that Ui is continuous and {ui > 0} H dU is con- 
nected. If 

^^^<Ai(l^) (13) 
s 

where Ai denotes the first eigenvalue of the Dirichlet operator with zero boundary condition, 
then {m > 0} is connected. 
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Proof: assume not, then there is a connected component A of {ui > 0} that does not 
touch the boundary; on the other hand, Uj > satisfies the equation 

-Aui = Ui 

\ Ui J 

on A with vanishing boundary trace. When testing the equation with Uj, by H13|) we obtain 
a contradiction with the Poincare inequahty. ■ 

It is worthwhile noticing that the condition H13j) is always satisfied on small domains. This 
can be useful in the local analysis of the solutions. Another useful property of S is that 
its elements are uniformly bounded in the interior of O thanks to the next proposition. 

Proposition 6.3 Let (ui) be an element ofS. Then the following hold 

(i) There are functions {^i) G W^^^i^) such that 

-A#i = fi{x,<^i) inn 

= Ui on on (14) 

$j > Ui in n, 

(a) There are functions (^'j) G W^^^{n) such that 

— A^'j = f{x,Ui) in 0, 

^i = Ui on do, (15) 

^i < Ui in 0, 

(in particular, "^f ^uf = Ui). 



Proof: to prove the first assertion we apply the method of upper-lower solutions: we 
need an ordered pair of functions < (3i where is subsolution and (3i supersolution of 
problem ((11)) . We simply let = Ui as lower solution; on the other hand we obtain a 
suitable /3j by solving 

-AA = hi{x)l3i in 

fii = M + (j)i on dVl 

for large constants M. Notice that assumption (A2) implies the existence of arbitrarily 
large positive functions /3j satisfying the above problem. Furthermore, since h{x)f3i > 
fi(x,Pi), then the (3i are supersolutions to equation (fTl|) . Finally we get (3i > Oi by the 
maximum principle. 

The proof of the second assertion is trivial, since the boundedness of Ui and the assumption 
(Al) imply the existence of a solution for problem (|15|) : the relation < Ui then follows 
by the maximum principle. 

To conclude we observe that the regularity of the <^j's and ^j's follows by the standard 
regularity theory for elliptic equations and our assumptions on the boundary data and the 
nonlinearities. ■ 
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Remark 6.2 As a consequence of the above proposition, the components of each element 
U S are uniformly bounded on compact subsets of 0,. Then, recalling (Al),(A2), if 
w CC r2 there exists M > (depending only on uj) such that 

U eS =^ -Aui < M, -Aui > -M on to (16) 

Furthermore, the regularity can be improved up to the boundary of in the sense that 
to bounded boundary data there correspond bounded barriers $i and ^j. Moreover the 
barriers will be Lipschitz continuous up to the boundary when both the data and the 
boundary di} enjoy the same regularity. 



7 The class S* 

Let M > and h he a fixed integer. We introduce 

Ui > 0, Ui ■ Uj = if i j I 
-Aui < M, -Aui > -M j 

It follows from Remark 16.21 that S C 5|f ^. It will be more convenient to work in this 
larger class rather that in the class S, for it is closed with respect to the limits of sequences 
of scaled functions. This property will be extremely useful in performing the blow-up 
analysis in Section |H1 We first present a technical result concerning with the elements of 
iS|^^(i?(0, 1)) when is a point of multiplicity at least two. 

Proposition 7.1 There exists M* such that, for all < M < M* the following holds 

(a) Let Un G Sl^j^{B(Q, 1)) such that m([/„)(0) > 2 and Un ^ U with U^O. Then 

Ui = l,...,h : Ui^0}>2. 



Definition 7.1 



(b) For every 7 > there exists C^j) > such that, if U £ S^.f i^{B (0,1)) with 
m(C/)(0) > 2 and |lf/||L2(^(-o^i)) > 7, then 

f [vt/p > C{j) f \U\^ 

JB{0,1) -^-8(0,1) 
where C > only depends on M* and 7. 

(c) LetUn G cSXj ;,(-B(0, 1)) such that m{Un){0) > 2; if \\VUn\\ 1^(3(0,1)) = 1 and \\VUn\\ 1^(33(0,1)) 
is bounded, then there exists U ^ such that Un ^ U. 
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Proof: Part (a). Arguing by contradiction we can assume the existence of M„ and 
Un e <S]^^_^(5(0, 1)) such that m(C/n)(0) > 2, Un ^ U with U ^ and N := ^{i = 
1, . . . ,h : Ui ^ 0} < 1. Note that, since U ^ 0, then must be equal to 1. Thus we 
can assume that U is of the form U = (ui, 0, . . . , 0) where ui ^ 0. Now consider (j) > 
solution of the problem —A(j) = 1 on -6(0, 1) with boundary conditions (p = 0. It holds 

-A(mi,„ + Mn(t>) > 

and by the mean value property for superharmonic functions we have 



= 2i,„,(0) > / Ui n{x)dx + Mn\ (l){x)dx - (/)(0) . (17) 

JBiO,!) ' \Jb(0,1) J 

By the compact embedding of in and since ui = ui, it holds 



ui n{x)dx — > / ui{x)dx = a > 0. 

5(0,1) ' "'-8(0,1) 

On the other hand, the last term of (|17)) becomes less than a/2 when M„ 0: this finally 
provides the contradiction > a/2. 

Part (h). Let 7 > be fixed. Assume by contradiction the existence of — and 
Un G ;,(5(0,1)) such that m(C/„)(0) > 2, \\Un\\L^i^B{Q,i)) > 7 and 

/b(0,1) |V^nP _^ Q 

/b(o,i) 

Let us define K := Un ■ ll^n||^2(B(o,i)): note that K G 5^^^^/^_^(S(0, 1)), m(T4)(0) > 2; 
furthermore ||V^||l2(5(o,i)) = 1 and || VV^||£,2(5(o,i)) ~* by construction. Then there ex- 
ists V such that Vn ^ V; moreover V ^ 0, since by the compact embedding of in 
it holds ||^||l2(-5(o,i)) = 1- Now note that, since the gradients vanish in the norm, then 
V = (ci, . . . ,Cfi) where Cj € M. Furthermore, since the supports of the components of Vn 
are mutually disjoint, passing to the limit a.e. in -6(0, 1), we get q • Cj = if i 7^ j: this 
means that only one of the components of V is not identically 0. This is in contradiction 
with Part (a) when applied to the sequence Vn- 

Part (c). By assumption \\^Un\\L^(^B{o,i)) is bounded: then, if M is small enough we can 
apply Part (b) and thus the whole -ff^-norm ||C^||_f/i(s(o,i)) is bounded. This provides U 
such that Un U. Now test the variational inequality — Anj_„ < M with Ui^n- 

[ |VUi „P < M [ Ui n + / Ui n-^Ui^n- 

JB{0,1) Jb{0,1) JdB{0,l) Of 

Assume by contradiction that U = Q: by the compact embedding of in and since 
II VC/„ 11^2(95(0,1)) is bounded, we deduce that the r.h.s vanishes. This implies ||VC/n||L2(5(o^i)) 
0, in contradiction with the assumption ||VJ7„|[j;^2(5(o^i)) = 1. ■ 

Remark 7.1 The argument used in the proof of part (a) allows to establish a mean value 
property for functions which are superharmonic up to a small term. Precisely, if (v„) C 
is such that —/S.Vn > —Mn with (0 <)Mn and /g(o i)Vn > a > 0, then Vn{0) > a/2 if 
n is large enough. 
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8 Lipschitz Regularity 



A key tool in studying the regularity of both the function U and the free boundary is a 
suitable version of the celebrated monotonicity theorem, see In this paper we shall 

take advantage of the following formula which is proven in [Sj . 

Lemma 8.1 Let wi,i«2 £ n such that —Awi < 1, • u;2 = a.e. and xq G 
d{{wi > 0}), i = 1,2. Then there exists C > 0, independent of xq, such that 



2 1 

IT l2 



\ywi{x)\'^ 

\N-2 



dx < C. 



JB{xo,r) \X - Xo\ 

In the proof of our regularity results we shall also need the following technical lemma. 
Lemma 8.2 Let U E H^{^) and let us define 



(18) 



1 



\\7U{y)\^dy. 



r'' JB{x,r)nn 
Lf {xn,rn) is a sequence in 0, x M+ such that 

(j){Xn,rn) OO 

then r„ — > and 

(i) there exists a sequence {r'^)n C such that (f){xn-,r'^) oo, and 



(19) 



,9 N 

V?7P < — 



aB{x„,r'jnn 



IVC/I 



r'n JB{x„,r'Jnn 



(ii) if A an and 



dist{xni A) 



< C 



then there exists a sequence {x'^,r'^) such that cj}{x'^,r'^) oo and x'^ G A for every 



n. 



Proof: since U G Ll^i^ and cj){xn,rn) +oo, obviously — 0. We begin proving (i). 
Let g he defined on the whole as 



g{x) :-- 



iv[/(x)|2 xen 
xeR^\n. 



Clearly 

{x,r)r\Q \'^U\'^ — JB[x,r) 9) IdB{x,r)nQ I VC/p — jQB{x,r) 9- We observe that 

N 



r^i \JdB{Xn,rn) 



9{y)dy 



Tfi JB{xn,r„) 



g{y)dy . 
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As a consequence, our problem is reduced to find such that (/)r(x„,r^) < and 
(/)(x„,r^) +00. Let := inf{r > r„ : (/)r(x„,r) < 0}. We have that < +00 
for every n (recall that (j){xn,r) ^ < (/)(a;„,r„) as r ^ 00), and (j){xn,rn) < (j){xn,r'^), 
and (i) is proved. In order to prove (ii), let G ^ such that dist{xn,x'j^) < 2Cr„ (this is 
possible by assumption) and let := (2C + l)r„. By construction B{x'^,r'^) D B{xn,rn), 
that implies (j){x'^,r'^) > (2C + 1)~^0(3;„, r„), and also (ii) follows. ■ 



8.1 Local Lipschitz continuity in S* 

The class S* seems to be the natural framework for proving the interior Lipschitz regularity. 
Indeed we have: 

Theorem 8.1 Let M > and k be a fixed integer. Let U G 5|^^(r2); then U is Lipschitz 
continuous in the interior offl. 

Proof: we consider the set ii' CC compactly enclosed in and the function (j){x,r) 
defined as in restricted to the set D := {x € n',r G M+ : 2r < dist{dn' ,dn)}. We 
have to prove that (p is bounded on D. We argue by contradiction, assuming that 

sup (/)(x, r) = +00 
D 

i.e. there exists a sequence (a;„,r„), such that 

rl^^^L miy)\'dy = +oo. (20) 

By Lemma |HI21(i), there exists a sequence (denoted again by r„) satisfying (|2U|) and more- 
over 

/ |VC/|2 < — f [VC/p. (21) 

JdB{ ) fn JB{xn,r„) 

Now we distinguish three cases, according to the nature of Xn and r^, up to suitable 
subsequences: 

CASE I: m{xn) = for all n (up to a subsequence) and > i 

We immediately obtain a contradiction with H20|) . since in this case Uj = on B{xn,rn), 
for all i. As a consequence, using Lemma [8.21 (ii). without loss of generality, we can assume 
that m{xn) > 1 for every n (besides (PT?]) and ((23)- 

CASE II: m(xn) = 1 for all n and > 1. 

In this case we can assume that only ui ^ on B{xn,rn) and we define 

Vn{x) = {Ui{x) - Ui(x„))+ Wn{x) = {ui{x) - Ui{Xn))~ 

Then w„ and Wn are in L°°{^}') by Remark 16. H they have disjoint supports, —Avn < M, 
— Awn < M and v„(a;n) = Wn{x) = 0. Thus we can apply the monotonicity formula of 
Lemma 18. II 

1 f \^Vn{x)\^ 1 f \VWn{x)\^ 



fn JB(xn,rn) \x Xn\^ ^ r"^ J B{xn,r„) \x ^ 
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uniformly. Since |x — x„| < r„, we deduce 



/ \VVn{x)f-—f \VWn{x)\^<C 

J B(xn.rT,) J B(xn.r„) 



where C is independent of n. Comparing with (|2U|) we have that only one of the two term 
is unbounded and forces the second one to vanish, e.g. 

1 /• ,„ ,0 1 



TV / I '^1 A' 

f-n JB{Xn,rn) ffl JB(Xn,rn) 



V7;„|2^oo ^/ \Vwn\'^Q (22) 



as n — > oo. Let 



2 1 /■ 

^" L, , 1^^" 

'n J B{Xn,rn) 



2 



and let us perform the blow up analysis around Xn with parameter L„ by defining Vn = 

{vi,n-,V2,n) aS 

1 1 

Vl,n{x) = J Vn{Xn + rnX) V2,n{x) = j Wn{Xn + TnX) XGB(0, 1). 

It is easy to verify that Vn G 5X/„ 2(-^(*-'' where M„ = rnM/Ln — > 0. By construction 
we have that /^(o i) l^^nP bounded; using (|21jl . this implies that JgB(oi) I^^^P is 
bounded too: thus Vn satisfies the assumptions of Proposition 17.11 This provides the 
existence of a weak limit V = (^1,^2) (by Part (c)) such that Vi ^ 0, i = 1,2 (by Part 
(a)). But this is in contradiction with H22() that forces f 2 = (by Part (b)). Again, this 
contradiction and Lemma 18.21 (11) allow us to assume m[xn) > 2. 

CASE IIL m[xn) > 2 for all n. 

In this case the proof follows closely the line of the previous one. Let us give some details: 
since m{xn) > 2, we can apply the monotonicity formula to each pair wi := Ui, W2 ■= uj.: 

1 f \Vwi{x)\'^ 1 f \Vw2{x)\'^ 



r2 / 7—) / \ I /y* /y" I W 2 'Y* 2 / T-j / \ I It* fy* I -/\^ 2 

n JB(Xn,r„) \X Xn\ 'n J B(Xn,rn) \X Xn\ 

uniformly. By (pn]) we deduce the existence of one index i such that, up to a subsequence, 
it holds 

-L/ iVn^i^^oo -1/ iVn.p^O VjV^ (23) 

'n JB{Xn,rn) 'n •'B{Xn,r„) 



as n — > 00. Let 



fn JB{x„,r„) 

and perform the blow up analysis around x„ with parameter L„ by defining T^i = (^^j,n) 
as ^ 

■"^•,^(2:) = J Uj{Xn + TnX) X G B{0, 1). 

It is easy to verify that Vn G 5^^^ ;j,(i?(0, 1)), where, again, Mn = rnM/Ln 0. By 
construction we have that /^(g i) l^^nP is bounded and, again by Jqb(oi) l^^nP 

is bounded too: thus Vn satisfies the assumptions of Proposition 17.11 This provides the 
existence of a weak limit V (by Part (c)) such that at least two of its components are 
strictly positive (by Part (a)). But this is in contradiction with ()23() that forces (by Part 
(b)) Vj = for all j ^ i. ■ 
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8.2 Lipschitz continuity up to the boundary in S 



In this Section we are concerned with the regularity of the elements of S up to the bound- 
ary, in the case of regular boundary and Lipschitz boundary data. Our main goal is the 
following result. 

Theorem 8.2 Let dQ be of class , U e S with Ui\Qn = (f>i and G Ty^'°°((90) for 
every i. ThenU eW^'^^iTL). 

The proof relies upon the local analysis as developed in the previous section, joint with 
the suitable use of the pinching property stated in Proposition 16.31 We begin with some 
preliminary remarks. 

Remark 8.1 Under the assumptions of Theorem XH.'A we have 

(i) Ui ^ C(r2) for every i (and, in particular, it makes sense to consider pointwise values 
of Ui); 

(n) ^laneL'^idQ) for every i. 

Proof: since U G S, through Proposition 16.31 we obtain the existence of /c-uples of func- 
tions {^i), (^i)> with the properties introduced in that proposition. Moreover, since 
G W^''^{dQ) for every i, by standard regularity theory for elliptic equations we infer 
(J). £ W^'°°{Tl), G W^'°^(U) for every i. By Theorem lO (which holds in this case 
through Remark Ui G C(r2); since "iff < Ui < and = on 50, (i) easily fol- 
lows. Moreover, using the same inequality and the very definition of directional derivative, 
we obtain, in distributional sense, 

d^t ^ duj ^ d^i 
dv ~ dv ~ du ^ 

and also (ii) follows. ■ 
Now we are ready to prove Theorem 18.21 

Proof of Theorem 18.21 let U satisfy the assumptions of the theorem. By the first part 
of the previous remark ui G L°°(r2) for every i, and hence there exists a constant M such 
that 

-Aui < f{x,Ui) < M. 

Moreover, by Proposition 16.31 there exist A;~uples (^i), such that "^f{x) < Ui{x) < 
^i{x) on $7, ^f{x) = Ui(x) = $i(x) on dQ and, as we just observed, ^i, G W^''^{^) 
for every i. 

As usually we define (pix, r) as in ()19() and we assume by contradiction the existence 
of a sequence {xn,rn) such that (j){xn,rn) +oo. Assume that (up to a subsequence) 
dtst{x„,dQ) y This means that B{xn,rn) H dO, = for every n. In this situation, one 
can repeat exactly the same proof of Theorem 18. II frouehlv speaking, in such a situation 
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the blow-up procedure does not "see" the boundary), obtaining the same contradictions. 
Hence we can assume that '^^■^*(^"'^^) < i a,]^^^ i^y Lemma l8?2l fii). without loss of generality 
we can choose Xn € dQ; moreover, we take r„ such that the inequality in Lemma l8.21 (i) 
holds. Let (/>(x„,r„) =: — +00 and define, for every i. 



1 



1 

1 



-{Ui{Xn + r„x) - Ui{Xn))., 

I 

-(^>,(x„ + r„x) - ^i{Xn)). 



We have that < Ui^^i^) < ^i,n{x) on 17^, where Q.n '■= {x € 5(0, 1) : x^+rnX G fi} 

(this inequality holds because the non-scaled functions coincide in x„ G dVL). We observe 
that, taking into account that r„ — > and L„ — > 00, we have 

-Anj,,i(x) = --^Anj(x„ + r„x) < -^-M < 1 
for every i, when n is sufficiently large. Testing the above inequality on Ui^n we obtain 



< 



dui 



— < ||'Wi,n||L°°{n„) + ll^j,ri||L2(9n„) 



dui 



dv 



Our aim is to prove that the righthand side of the previous inequality tends to for every 
i. This will provide a contradiction with the fact that, by construction, J2 In„ = 1 

for every re. 

Clearly ^'i,n(0) = $i,„(0) = 0. Moreover, 



|V$i,„||LOO(Q^) 



Lr 



\^^i\\L°°{nnB{x„,rn)) 0, 



and the same holds for ^i.^. This implies ||*J>j,n|lvyi.°o(o„) ^ 0^ ~^ ^ 

therefore 

ll«j,n||L°°(n„) ^ Vi 

Since dO, is of class C^, we obtain that d^n has bounded (A^ — 1) -dimensional measure, 
and thus also 



Wi,n\\L^{dn„) Vi. 



I dui „ II 



is bounded. To this 



Therefore the only thing that remains to prove is that 
aim, let d^n = Ti.n U T2^n, where Li^^ := dB{0, 1) n r2„ and T2^n '■= dVtn \ ^i,n- Then the 
estimate on Fi^^ descends from Lemma l8.2( (i). recalling that it implies 









dv 



< 



r„ JdBix„,r„)nn 



N 



N. 



On the other hand, the estimate on T2^n is an easy consequence of the bounded measure 
of r2,n and of E,emark IS.ll fiil. ■ 
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9 Further regularity in dimension N=2 



9.1 Vanishing of the gradient at multiple points 

Let N = 2 and U G S. The main goal of this section is to prove that the gradient of U 
vanishes continuously at points of multiplicity at least three. This result will be established 
through the application of a monotonicity formula with three or more phases. To start 
with, we need the following technical result, that allows us to reduce the Ui^s to solutions 
to suitable divergence-type equations. 

Lemma 9.1 Let a € and let v be an solution of —Av < a{x)v in B{xo,r). 
Then, if r is small enough, there exists (p G such that (p is radial with respect to xq, 
infs(xo,r) > and 

-div((^2(V-)) < in B{xq, r). 

Proof: let us consider the eigenvalue problem 

— Au(x) = a(x)u(x) x G B{xo,r) 
u{x) > X G B{xo,r) 

If r is small enough, the above problem can be solved in the class of radial functions with 
respect to xq. Let (p be such a solution: now by elementary computations 



-if^A 2ipV-Vip < 

ip ip 



giving the required inequality. 



This local reduction will be widely exploited throughout the present and the next section. 
As a first application it allows to prove a variant of the original monotonicity formula by 
Alt-Caffarelli-Friedman . 

Lemma 9.2 Let U G S and Wi = J2jeii ^j; where Ii L) ... U Lh C {1, k}. Assume that 
xq G d{{wi > 0}). Then for all h >2 there exists C , independent of xq such that 

\\\ I \Vwi{x)\''dx<C. (24) 

f^^ r^ JB{xo,r) 

Proof: (sketch). Let xq G d{ui > 0} and recall that, since [/ G 5, it holds —Aui < fi{ui) 
for all i. Let 0(2;) := maxj=i ^ note that a G L°° by the assumption on fj. 

By Lemma l9. 11 there exists r > and a regular radial function ip which is strictly positive 
on B{xo,r) such that, for all i 

-dW{ip'^(y—)) < in B{xo,r). 



21 



Now consider Wi = ^ji^i^Uj, where Ii U ... U I/i C and let Wi := Wi/ip: then 

-(MY{<f'^{Vwi)) < 0. Let = 2 and set 

^ I f 

= n ^ / ^'^{x)\Vm{x)\^dx. 

Following the proof of Lemma 5.1 in |2] we can compute 

$'(r) > f-Zi+^infV yX;"^ 

where P = {(ri)f^i, UL^ = 55(0, 1), meas(ri n Fj) = 0} and 

Jr. \dev? 



A,- = inf 



2 ■ 



Since Aj > (j^) we immediately obtain that <I>'(r) > in [0,r']. 

By this formula and since there exist positive a < b such that a < ip{x) < b for all 
X G B{x(),r'), it follows in particular that 1)24^ holds. ■ 

Now we are ready to prove the main result of this section. 



Theorem 9.1 If xq E Z-^ then \VU{x)\ — > as x — > xq. 



Proof: assume by contradiction the existence of Xn — > xq and r„ ^ such that 

lim f |VC/p = a (25) 

a G (0, oo]. 

We first claim that equation (|25)) holds for Xn G ^3- The proof of this fact can be done 
as follows: let p, r > be fixed and let Ap^r = {x G B{p,xo) : d{x, Z^^U)) > r}. By 
Proposition 16. H in Ap^r we can give alternate positive and negative sign to the Uj's in such 
a way that the resulting function v locally solves an equation of the form —Av = f{x, v) 
where f{x,v{x)) = fi{x,v{x)) if v{x) = ±Ui{x). For x G Ap^r we define 

$(x) = 4/ \VU{y)\^dy = \ f \VU{x + y)fdy. 

JB{x,r) JB{0,r) 

By elementary computations — A(|Vf (x)p) = —2\Vv{x)\'^—2fx{x, v)\Vv{x)\—2fs{x, f )))|Vf (x)p: 
recalling that v is bounded by Remark 16.21 and since |VC/[ = \'Vv\, this gives — A$ < a$ 
on Ap^r, for some positive constant a (depending on r). Then, by an extension of the max- 
imum principle, there exists C > (independent of r) such that max^^ ^ $ < C maxg^^^ 
This implies that for n large enough, H25() holds with a = a/2 and for a choice of such 
that d{x'n, Z^^U)) X r„; then taking z„ G Z^{U) such that |z„ — x^j = d{x'^,Z^{U)) we 
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obtain that l|25() holds for balls centered at z„ and radius r„ + d{x'n, Z3{U)) x r„. 
Furthermore, by exploiting the Lipschitz regularity of U, we have 

/ lV[/p < — / [VC/p (26) 

JdB(x„,Pn) Pn J B(x„,pn) 

where pn = {2L/C)rn and L is the Lipschitz constant in 17' = {x € : d{x,dil.) > 
d{xo,dn)/2} D {Xn}. 

Hence, in the following let us assume that 1)251) holds for a choice of Xn G 2^ and of radii 
Pn satisfying (we denote pn again by r„). 

Let r > 0: by the monotonicity formula with three phases, there exists C > (independent 
of r) such that 



3 

f] \ f \Vuj,{x)\'^dx < C 

fj^ JB(x„,r) 



for all LOi := Ui, L02 '■= Uj and los := J2h^{i,j} ''^h, such that x„ € (?{nj > 0} n d{uj > 
0} n d{ui > 0} for some / {^j j}- Then, due to assumption (|25|). we deduce that there 
exist at most two components, say ui and U2, such that 

lim ^ / |VMi(x)p(ix > 0, i = 1,2 4-/ |Vui(2;)pdx ^ V i > 3. 

(27) 

Let us set 

Ll:=\ [ \VU{x)\^dx 



JB{Xn,rn) 

and consider the sequence of functions 

Unix) = -^U{xn + r„x) (28) 

defined in x G 5(0,1). Note that, since U G S^j by (Unj, then [/„ G ^ where 
Mn := f^M. Then /^(^ -^^ |VC/„p = 1 and, by JgB(o,i) |Vf/n.|^ is bounded too. By 
Proposition 17.11 Part (c), there exists U G H^{B{0,1)) such that, up to subsequences, 
Un ^ U and C/ 7^ 0. Furthermore, by Part (a), we know that at least two components of 
Un does not vanish to the limit: comparing with ()27|) we have 

U = (ui,U2,0, ...,0) 

with Ui ^ Q for z = 1,2. Moreover, since M„ 0, it holds U & Sq^. This in particular 
implies that ui — U2 is harmonic. 

Now, if ui(0) > 0, (resp. 112(0) > 0) then ui > (resp. U2 > 0) in B{0,f) for some f > 0. 
This implies Jb{o f)^^,n{x) > a for some a > 0, since it converges to the L-'^-norm of lii 
in B{0,f). We can thus apply Remark l7. II to the sequence (ui^n) obtaining ui,n(0) > a/2. 
But this is in contradiction with the fact that, by definition, ui^niO) = 0. 

Now set V = ui and assume that v{0) = 0. Standard results on harmonic function (see 
[21\ ) imply that v{r, 9) ~ cosp{6 + ^o) for some p > 1. Thus, by the strong convergence 
in and a diagonal process, we can assume that, for n large enough and i = 1,2, there 
exists rrii > such that Un,i > mj on a circular sector An,i = {{p,9) : r < p < R, an,i < < 
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Pn,i} C {un,i > 0}. Here we can assume, for instance, an,i = a < Pn,! < Q<n,2 < Pn,2 = P 
with a and /? fixed; note that a„^2 — ^ as n ^ cxd. Now, since is a zero of Un with 
multiplicity m(0) > 3, there exists a third component, say n„ 3, and a continuous path 
7n : [0,1] ^ S(0,i?) such that 7„(0) = 0, 7„(1) G dB{0,R), 7„(t) G < < a„,2}, 

'"n,3(7n(*)) > for all t e (0, 1]. Therefore, if we set S = {{p,9) : r < p < R,a < < (3} 
and denote by LOn,i {i = 1,2) the connected component of {un,i > 0} that contains An^i, 
then dist{S fl iOn,i, S n u;n,2) > 0. 
Now consider 

and for x G Tn{S) define 



We claim that 



Ui{x), xeTn{Un,i) 1=1,2 

-Ui{x), x G {uj > 0} \ (r„(a;n,i) ur„(a;„,2)) i = l,...,k 

- > /K) in Tn{S) (29) 



where f{x, s) := fi{x, s) for x G {ui > 0}. 

In order to prove this assertion, let us drop the dependence on n. Then fix any (j) > 0, 
(f) G Cq{T{S)): we have to prove 



X 



{VwVct)- f{w)(t))dx>Q. 

T{S) 

For easier notation, set Ai = T{uji), A2 = T{u2) and B = T{S) \ {Ax U A2). Now take 
a partition of the unity in such a way that (l){x) = 4>i{x) + (f>2{x) where (f)i > and 
> 0} C AiU B. For x G r(S') define the function 



wi{x) 



ui{x), X e All) {B n{ui > 0}) 
—Ui{x), X E {ui > 0} i = 2,...,k 



(Analogous definition for W2-) Since U E S (namely — Aui > /(ui)), it holds 

(VwiV^i - /(wi)0i) dx>0 



I 

Jt{ 



IT{S) 

and 



/T(S) 

Summing up the two inequalities we obtain 



/ {Vw2V(j)2 - f{w2)(t>2) dx > 0. 

Jt(s) 



f (VwVcl) - f{w)4>) > -2 V / {VuiV(l)i - f{ui)(l)i) dx. 

Jt{S) -t^JB 
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As we are going to prove, each term of the above sum is negative, and this finally proves 
assertion ()29p. To this aim let i = 1, and introduce a cut-off function rj such that rj = 1 
on 5 n {til > 0} and ?? = on T{S) \ B. Let tp = ri4>i e C^{T{S)); then 

/ (VniV0i-/(ui)0i)dx= / {VuiViP- f{ui)iP)dx. 
Jb Jt{s) 

Then the righthand side is negative since U & S (namely —Aui < f{ui)). 

Final step. Now fix t„ > and > r„ > r in such a way that, if we set y„ = ^n{tn) then 
it holds jQB{yn r„) > « for some positive a. (This is due to the fact that dB{yn,rn) C 
LOn,! ^i^n,2 except for a small piece of total length less then 2i?(an,2 — /3n,i) — J- as n ^ oo.) 
This allows to apply Remark l7. II to the sequence of rescaled functions ( ^"'£ ) (as in (^5]) ). 
since they satisfy A ^"*]- > — M„ 0. This gives WniUn) > 0, in contradiction with the 
fact that, by construction, Wn{yn) = — ^n,3(yn) < 0. ■ 

9.2 Local properties of the free boundary 

Let again N = 2; the purpose of this section is to deepen the investigation of the behavior 
of the elements U of the class S around multiple points, that is, x € with ■m{x) > 2. 
Our main goal is to prove that, near a multiple point, U and its null set exhibit the same 
qualitative behavior of harmonic functions and their nodal sets. We refer the readers to 
the fundamental papers of Alessandrini 1 and Hartman Winter j21j for the main results 
about the zero set of harmonic functions and, more in general, of functions in the kernel 
of a divergence type operator. We shall obtained the desired description for those U (z S 
with the property that each component has connected support, i.e. 

U G S such that {ui > 0} is connected (30) 

In order to simplify the topological arguments involved in the proof of our results, in what 
follows we shall make the additional (but not necessary) assumptions that is simply 
connected and that all the boundary data are not identically zero, namely ^ on dil.. 
In this situation the assumption ()3U() (and Remark 16.1(1 implies that {ui > 0} is simply 
connected. 

Let us start our description by analyzing points of multiplicity two, that will be denoted 

by 

= {x eft: m{x) = 2}. 

We already noticed in Remark 16.11 that this set is locally a regular arc around those 
points where VU does not vanish. Our first result states that in fact this is always true: 

Lemma 9.3 Let xq ^ Q such that m{xo) = 2. Then VU{xo) ^ and Z'^ is locally a 
-curve through xq. 

Proof: since xq G d{ui > 0} fl d{u2 > 0}, then for all r small enough B{xo,r) n {ui > 
0} = for all i > 2. Assume by contradiction Vui{xq) = = Vu2{xo); then u = ui — U2 
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is and satisfies the equation —An = a{x)u (with a{x) = ^^('j.jZu^^") ) in B{xo,r). 
By Lemma l9. II there exists a positive, regular function (p which is radial with respect to 
xq such that 



on B{xo,r). Then u/(p satisfies all the assumptions necessary to apply the main theorem 
in which says that the null level set of u/ip (indeed of u) near xq is made up by a 
finite number of curves starting from xq. Obviously in our situation such number must be 
even. Now recall that each {ui > 0} is connected in fi: by a geometrical argument we can 
see that the null level set is made up by (two semi-curves joining in) one C^-curve. But 
again applying [T] we have VU{xq) 7^ 0, a contradiction. ■ 

Lemma 9.4 Let xq Z3. Then there exists {xn} C il, such that m{xn) = 2 and x„ — ^ xq 

Proof: assume not, then there would be an element uq of Z3 having a positive distance d 
from Z'^. Let r < d/2: then the ball B{yo,r) intersects at least three supports; therefore 
there exist, say, x G {ui > 0} and zq G Z3 such that p = d{x,zo) = d{x,Z^) < d{x,Z'^). 
Then the ball B{x,p) is tangent from the interior of {ui > 0} to Z^ in zq; furthermore Ui 
solves an elliptic PDE and it is positive on its support: we thus infer from the Boundary 
Point Lemma that Vui{zQ) 7^ 0, in contrast with Theorem 19. II ■ 

Let us now prove an asymptotic formula describing the behavior of J2 ™ the neighbor- 
hood of a multiple point which is isolated in Z3. 

Theorem 9.2 Let xq € Q with m(xo) = h > 3 that is isolated in Z3. Then there exists 
9q G (— vt, tt] such that 



as r ^ 0, where {r,9) denotes a system of polar coordinates around xq. 

Proof: by assumption xq is isolated in Z3, then there is B = B{xq, p) such that Z^CiB = 
{xq}; furthermore since each {uj > 0} is simply connected, then {uj > 0} n dB = 6j is a 
connected arc on dB for the h indices involved in xq. Choosing a slightly smaller radius we 
can suppose that the intersection of dB with Z'^ is transversal (recall that, by Lemma 19.31 
Z'^ is locally a C^-curve). Let us assume that h is even: then we define a function v{r, 9) 
such that \v{r,9)\ = \u{r,9)\ and sign{v{r, 9)) = ( — I)-' if (p, ^) (zbj,j = l,...,h. Note that 
the resulting function is alternately positive and negative on the consecutive (with respect 
to 9) local components of U. If on the contrary h is odd, we define |f(r, ^)| = \u{r'^,26)\ 
and we prescribe an alternating sign to the local components of u{p'^, 26). It is worthwhile 
noticing that the resulting function v is of class in i? = B{xq, p^): indeed, {ui > 0}ni? is 
simply connected for every i and thus each connected component of B\u~^{{)) corresponds 
to two components of i? \ ^"^(0) to which we give opposite sign. In both the even and 
the odd cases v is of class and it solves an equation of type — At> = a{x)v in B\ {xq} 
(resp. i? \ {xq}), where a € L°° is given by f{v)/v and r'^f{v)/v respectively. Moreover 




U{r,9) =rt|cos(-(6' + 6'o))| + o(r2) 



26 



Vw(xo) = by Theorem 19.11 this implies that v is in fact solution of the equation on 
the whole of B (resp. B) and thus it is of class C^'°. Now, choosing r small enough, by 
Lemma l9. II we have a positive, regular function ip such that 

-d„(,^v(^))=o. 

Then, we can apply to v/(p the asymptotic formula of Hartman and Winter as recalled 
in . To complete the proof, let us observe that h represents the number of connected 
components of Li{ui > 0} in a ball centered in xq, providing the choice of the periodicity 
of the cosine in the representation formula. ■ 

The last part of this section is devoted to prove that consists of (a finite number of) 
isolated points: thus Theorem 19 . 21 will provide the complete description of U near multiple 
points. 

We start by proving an intermediate result: 

Proposition 9.1 has a finite number of connected components. 

Proof: let := {uj > 0}; take an index pair (i, j) such that duji, duj do intersect and we 
consider 

Tij = duJi n dujj n z'^ 
uJij = uJiU ujj U Fjj- . 

Since by Lemma 19.31 r, j is locally a regular arc and each loi is open, it easily follows that 
Ui^j is open. Furthermore is simply connected. Indeed, let us consider a loop 7 in 
u)ij which is not contractible in Uij. This means that the internal region of the loop, 
contains at least a multiple intersection point. Since each support is connected, there 
exists w/j, h 7^ i,j such that lo]^ C ^2-^. But this is in contradiction with the fact that 
ojhnd^ = {(t)h > 0}. 

A similar reasoning allows to prove that each Ti^j consists in a single C^-arc: as straight- 
forward consequence of these facts the set of multiple points can have only a finite number 
of connected components. ■ 

Now we will need the following definition of adjacent supports: 
Definition 9.1 We say that uji and loj are adjacent if 

Let us list some basic properties: 

1. Every uJi is adjacent to some other ujj. This follows from the Boundary Point Lemma. 

2. Let us pick k points S uJi, i = 1, . . . ,k. If and ojj are adjacent, i < j, then 
there exists a smooth arc jij with 7ij(0) = Xi, 7ij(l) = Xj lying in coi Uuj U {vij}, 
for some i/ij € Z'^. 
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3. We can choose the arcs ^ij in a manner that they are mutually disjoint, except for 
the extreme points. 



We call Q the graph induced by the arcs ^ij and their endpoints. 

Construction of an auxiliary function v. Aim of this paragraph is to build up by the 
components oiU function carrying a sign law which is compatible with the adjacency 
relation and solves an elliptic equation. Let us first assume that 



In this situation the graph is a disjoint union of a finite number of branches: we select one 
of them. Now we define v as follows: if the index i is not involved in the branch, we set 
t; = in {ui > 0}. Next we follow the selected branch of the graph and prescribe a sign 
to each vertex by alternating plus and minus. We then define v{x) = ibii(x) according to 
this sign rule. Taking into account Remark 16.11 and Theorem 19.11 with this procedure we 
obtain a C^-function on 

Moreover v solves an elliptic equation as follows: let us define f{x,s) := fi{x,s) if a; G 
{ui > 0}, i = 1, . . . , k. Then by construction 



where a = f{v)/v G In fact we are going to prove that v solves the elliptic equation 
on the whole of 0,. We need a technical lemma 

Lemma 9.5 Let = {x Q : m(x) = 2} n Bi,[Z^). For every (5 > there exists e > 
such that 



Proof: by testing the equation (|31|) with the test function ip = 1 and integrating over the 
set Ui > a we obtain the bound, independent of a and i. 



Q has no loops. 



Av = a{x)v 



in Q \ 2^3 



(31) 





and therefore, passing to the limit as a ^ 0, 




The assertion then follows from Lemma l9.4| together with Proposition 19.11 



Lemma 9.6 v solves —Av = a{x)v in Vl. 
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Proof: let us fix a connected component of 2^3, named A. Thanks to Lemma 19.41 for any 
e we can take a neighborhood of A, say Ve C Bf;{Z^), in such a way that the boundary 
dVe is the union of a finite number of arcs of Z"^ and supplementary union of pieces of 
total length smaller than Ce. Let be a test function. We write 

j^{VvVLp- a{x)vLp) = 

= /q\v,(^^VV' - a{x)vLp) + /v_jVfV99 - a{x)vip) < 

<C^/avjVH + C/v.(|VH + l^l)- 

Let (5 > 0: we can find e > such that Lemma l9 . 51 holds . Moreover, by Theorem 19.11 we 
can assume that e is taken so small that supy^(|Vf | + \v\) < 5. Hence the above integral 
is bounded by C6. Since 6 was arbitrarily chosen we obtain that v solves the equation in 
a distributional sense. Usual regularity arguments allow us to complete the proof. ■ 

With this we easily deduce the final result, that is, U has only a finite number of multiple 
points: 

Lemma 9.7 The set Z^ consists of a finite number of points. 

Proof: since by Lemma 19.61 v is a solution of —Av = a{x)v in 0, locally we can reduce 
u to a function in the kernel of a divergence-type operator as in Lemma l9.ll Then the 
results of ^ and [U ensures that v has only a finite number of multiple points. ■ 

Let us now go back to the case when the graph associated to U presents loops: 

G has a loop. 

Let us first define an order relation between loops, according whether one is contained 
in the interior region of the other. Let us select a minimal loop 7 (no other loops are 
contained in Q^, its interior region). We can assume that contains at least an element 
xo € Z3 (if not, we can perform a conformal inversion exchanging the inner with the outer 
points): fix xq as the origin. Note that the supports involved in a minimal loop have the 
remarkable property that, for all i, oJi is adjacent to coj for only two indices different from 
i. Thanks to this property, if the number of vertex of 7 is even, we manage to assign a 
sign law to all the subset of Q contained in so that adjacent supports have opposite 
sign: it suffices to follow the loop and prescribe alternating sign to its vertices. We define 
v{x) = =bn(x) according to this law. 

If the number of vertex of 7 is odd, we wish to "double" the loop. To this aim, we define 
new Wj's by taking the complex square roots of the old ones. In this way the new loop 7 
will have an even number of edges and we define ?;(r, 9) = ibit(r^, 20) by giving alternating 
sign at the vertices of the loop. In this way we define a function v in $7^ which is of class 
thanks to Remark 16 . 1 1 and Theorem 19.11 Again, v solves the elliptic equation 

— All = a{x)v in 0^ \ ^3 (32) 

where a is either a = f{v)/v or a = r'^f{v)/v according to the construction of v. As in the 
previous case, it is possible to show that v solves the equation on the whole of fi^. The 
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proof follows exactly that of Lemma 19.61 with two remarks; first note that Lemma l9.5l still 
holds. Furthermore, in this situation the set Zj, n 0^ is connected (assume not; then in 
there are two points of Z-^ which are connected by a regular arc of double points. By a 
simple geometrical argument in the plane this implies that one of the supports is adjacent 
to other three different supports: as already observed, this is in contradiction with the 
minimality of the loop). Then will be a neighborhood of the whole ^3 H il^ with the 
properties required in the proof of that lemma. 

This immediately provides 

Lemma 9.8 There is only one point of Z^ lying in the interior of a minimal loop 7. 

Proof: let Vtr^ denote the internal region of 7: since by Lemma 19.61 u is a solution of 
— = a{x)v in ^7, locally we can reduce u to a function in the kernel of a divergence- 
type operator as in Lemma l9. 11 Then by and j21j we know that only a finite number of 
points of Z2, lie in Vt^. On the other hand, by the minimality of 7, Z^, H Vt-y is connected. 
Thus the origin is the only one multiple point contained \n.VL^. ■ 

We can finally prove that multiple points are isolated 
Theorem 9.3 The set Zj, consists of isolated points. 

Proof: recalling Proposition 19.11 we argue by induction over the number h of connected 
components of the set Z^. If /i = 1 then, by Lemma l9. 81 there is at most one minimal loop 
of the adjacency relation. If there is one, then Lemma 19.81 gives the desired assertion. If 
there are none, the thesis directly follows by Lemma [9. 71 Now, let the Theorem be true for 
h and assume that Z^ has h + 1 connected components. Again, if the adjacency relation 
has no loops we are done. Otherwise, we apply Lemma 19.81 to treat those connected 
components contained in the interior of the minimal loop and the inductive hypothesis to 
treat all those contained in the outer region. ■ 

Remark 9.1 Having proved that the multiple points are isolated, the existence of points 
of multiplicity zero can be easily ruled out for connected domains. 
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